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ABSTRACT 

Hard X-ray imaging of the Galactic plane by the INTEGRAL satellite is uncovering large numbers of 20- 
100 keV "IGR" sources. We present results from Chandra, INTEGRAL, optical, and IR observations of 4 IGR 
sources: 3 sources in the Norma region of the Galaxy (IGR 116195^945, IGR J 16207-5 129, and IGR J16167- 
4957) and one that is closer to the Galactic center (IGR J17195-4100). In all 4 cases, one relatively bright 
Chandra source is seen in the INTEGRAL error circle, and these are likely to be the soft X-ray counterparts of 
the IGR sources. They have hard 0.3-10 keV spectra with power-law photon indices of F = 0.5-1.1. While 
many previously studied IGR sources show high column densities (A^h ~ lO^^*"^"* cm~^), only IGR J 16 195-4945 
has a column density that could be as high as lO^^* cm~^. Using optical and IR sky survey catalogs and our 
own photometry, we have obtained identifications for all 4 sources. The /-band magnitudes are in the range 
14.9-10.4, and we have used the optical/IR spectral energy distributions (SEDs) to constrain the nature of the 
sources. Blackbody components with temperature lower limits of >9400 K for IGR 116195^945 and > 18,000 
K for IGR J16207-5129 indicate that these are very likely High-Mass X-ray Binaries (HMXBs). However, for 
IGR J 16 167-4957 and IGR 117195^100, low extinction and the SEDs indicate later spectral types for the 
putative companions, indicating that these are not HMXBs. 

Subject headings: stars: neutron — X-rays: stars — infrared: stars — stars: individual (IGR J16195^945, 
IGR J 1 6207-5 1 29, IGR J 1 6 1 67-4957, IGR J 1 7 1 95-4 1 00) 



\ 1. INTRODUCTION 

The hard X-ray imaging of the Galactic Plane by the INTEr- 
, national Gamma-Ra y Astrophysics Laboratory (INTEGRAL, 
■ lWinkleretalJl20()l ^ is uncovering a large number of new 
or previously poorly studied high energy sources. During 
the first two years of INTEGRAL operations (2002 October 
- 2004 September), 56 new "IGR" sources were discovered 
teird et al.ii 2006). and many more IGR sources have been dis- 
covered to date. INTEGRAL is also detecting many sources 
that are present in other X-ray catalogs but were not targets of 
focused studies until they were shown to be strong emitters of 
hard X-rays by INTEGRAL. If these sources are included, well 
over 100 IGR sources have been found*". All or nearly all of 
the IGR sources have been found in the 20-50 keV band with 
; the Imager on Board the INTEGRAL Satellite (IBIS) coded 
' aperture mask instrument iUbertini et al.ll2003t iLebrun et al.l 
1^03). Large numbers of IGR sources have been found be- 
cause of the combination of hard X-ray imaging with 12' an- 
gular resolution, a large field of view (9° by 9° fully coded 
FOV for IBIS), and the INTEGRAL observing plan, which 
emphasizes observations of the Galactic Plane. 

Follow-up observations of the IGR sources have shown a 
diversity of source types, including Low-Mass X-ray Bina- 
ries (LMXBs), High-Mass X-ray Binaries (HMXBs), Active 
Galactic Nuclei, Cataclysmic Variables as well as other source 
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types. While some of the sources have proved to be tran- 
sient, others have been consistently detected in X-rays in mul- 
tiple INTEGRAL observations as well as by other current and 
previous X-ray instruments, suggesting that they are persis- 
tent. One sub-group of IGR sources exhibits persistent but 
highly variable X-ray emission. These sources lie within a 
few degrees of the Galactic Plane an d show some evidence fo r 
clustering near Galactic spiral arms JLutovinov et al.ll2005at . 
Their X-ray fluxes are typically ~1-10 millicrab in the hard 
X-ray band (iDean et alJ 1200.5). and their luminosities are 
mostly unknown due to large uncertainties on their distances, 
but they may have luminosities of lO-'-'"'"' ergs s"' if they are 
at typical 1-10 kpc Galactic distances. 

Although X-ray, optical, and infrared (IR) observations of 
this group of persistent. Galactic IGR sources show that they 
do not have uniform properties, some trends have been iden- 
tified. For many of these sources, their X-ray spectra show 
high column densities with values of A^h well in excess of 
the levels expected due to interstellar material. The most ex- 
ti-eme example is IGR J16318-4848, for whichA^e 2 x lO^^ 
cm"2 (Matt &_Guainazzi 2003; Walte r et al. 20Q3), and there 
are also several ot her sources with values of A^h in the 10^^"^^ 
cm"^ range ( e.g., 'R odriguez etani2003t ICombi et al.ll2004t 
iBeckmann et al. 2005|). Many of the Galactic IGR sources 
also exhibit X-ray pulsations, indicating the presence of a 
neutron star Pulsations are detected for at least 5 of the 
persistent IGR sources with periods of 139- 1303 s (e.g., 
iLutovinov etlD l2005bt fBoda^hee et all l2006h . Finally, in 
some cases for which optical or IR spectra of IGR sources 
have been obtained, high-mass stellar companions have been 
found. IGR J 163 18-4848 harbors an unusual supergiant B[e] 
star, and the IR spectra also show P Cygni p rofiles, suggest- 
ing th e presence of a strong stellar outflow (iFilliatr e & Chatvl 
I2OOI . Another example is IGR 117391-3 021 (= XT E J1739- 
302), which also contains a supergiant iSmith et alJ l2006t 
^egueruela et al. 2006). 

Here, we present results for X-ray, optical, and IR follow- 
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up studies of 4 IGR sources that are in the Galactic Plane. 
Three of the sources (IGR J16195-4945, IGR J16207-5129, 
IGR J16 167-4957) are in the Norma region of the Galaxy 
(i.e., close to the tangent to the Norma spiral arm) at 332° < 
/ < 334°, and the fourth (IGR J17195-4100) is about half- 
way between the Norma region and the Galactic center at 
I = 347°. Although none of the sources have been previously 
well-studied, it has been pointed out that IGR J 16 195^945 
is probably the Advanced Satellite for Cosm ology and As- 
trophysics (ASCA) source AX J 16 1929-4945 ('Sugiza ki etaP 
UpOl; SidoH et al. 2005). The 2-10 keV ASCA spectrum of 
this source is well-described by a very hard power-law with 
a photon index of F = 0.6;!;q 5 and a relatively high column 
density of A^h = (12!^) x 10^^ cm . Based on its X-ray prop- 
erties, it has been su ggested that the source is a neutron star 
HMXB (ISidoU et aL 2005; Bird et al. 2006). IGR J16167- 
4957 and IGR J17 195^100 have been tentatively identified 
with ROSAT sources (IStephen et aljl200 5). but it appears that 
INTEGRAL is the first X-ray satellite to detect IGR J 16207- 
5129. The nature of these 3 sources is completely uncertain 
( iBirdet al. 2006) . 

In this wor k, we use observations with the Chandra X-ray 
Observatory dWeisskopf et alJl2002l) to obtain sub-arcsecond 
X-ray positions for the 4 sources (§2) and to determine 
whether the sources are associated with known sources (§3). 
In §4 and §5, we present a study of the Chandra and INTE- 
GRAL spectra. In §6, we report on optical and IR photometry 
of these sources we made at ESO's New Technology Tele- 
scope (NTT), and we combine these measurements with those 
available in optical and IR catalogs in order to constrain the 
optical and IR spectral energy distributions (§7). In §8 and 
§9, we discuss our results and present our conclusions. 

2. CHANDRA OBSERVATIONS AND SOURCE DETECTION 

We obtained short, ~5 ks, Chandra observations of the 
fields of IGR J 1 6 1 95^945 , IGR J 1 6207-5 1 29, IGR J 1 6 1 67- 
4957, and IGR J17195-4100 with the primary goal of local- 
izing the sources to facilitate IR and optical identifications. 
We chose targets located close to the Galactic Plane that were 
detected in the 20^0 keV band during our 2003 February IN- 
TEGRAL ob servation of the black hole transient 4U 1630^7 
iTomsick et al. 2004, 2005) and that were also rep orted in the 
first and sec ond catalogs of INTEGRAL sources (iBird et al.l 
I2OOI l2006l) . Consistent hard X -ray detections gave us rea- 
son to believe that these were persistent sources and that the 
probability of detection with Chandra was high. 

Table [2 shows the Observation IDs (ObsIDs) and times of 
our Chandra observations of the 4 IGR sources, which took 
place between 2005 April and Jul y. We used the Adv anced 
CCD Imaging Spectrometer (ACIS. lGarmire et al.'20'03) with 
the center of the 3' INTEGRAL error circle placed at the nom- 
inal aimpoint for the ACIS-I array. Although the effective 
exposure times are lower in the gaps between ACIS-I chips, 
we used a dithering pattern with an amplitude twice as large 
as the standard pattern to achieve a more uniform sensitiv- 
ity over the INTEGRAL error circle. In processing the data, 
we started with the "level 1" event list produced by the stan- 
dard data processing with ASCDS versions between 7.5.3 and 
7.6.2. For further processing, we used the Chandra Interactive 
Analysis of Observations (CIAO) version 3.3.0. 1 software and 
Calibration Data Base (CALDB) version 3.2.1. We used the 
CIAO routine acis_process_events to obtain a "level 
2" event list and image. 

The Chandra images are shown in Figure[2 For each target. 




Fig. 1. — Chandra 0.3-10 keV images of the 4 IGR sources. The large 
circles are the 3' INTEGRAL error circles. In three cases, ASCA or ROSAT 
counterparts have been suggested, and their error circles are also shown. Ar- 
rows point to the brightest Chandra source detected in each observation, and 
these are very likely the soft X-ray counterparts to the IGR sources. The 
lengths of the N and E arrows are l'. 



we searched for sources in a 7 x 7 arcmin^ region that includes 
the full INTEGRAL error circle. We restricted the energy 
range to 0.3-10 keV and used the CIAO routine wavdetect 
Uireeman et al. 2002). Based on the image size of 854-b v-854 
pixels and our detection threshold of lO"*" (see Freema n et alJ 
a precise definition), we would expect to detect < 1 
spurious source. For ObsIDs 5471, 5472, 5473, and 5474, 
we detected 4, 7, 4, and 7 sources, respectively. The bright- 
est sources for these four ObsIDs have 183, 678, 866, and 
676 counts, respectively, while the remaining sources have be- 
tween 3 and 24 counts. In addition, each of these four bright- 
est sources are located within the respective INTEGRAL error 
circles (see Figure[3, making it very likely that these sources 
are, in fact, the soft X-ray counterparts of the IGR sources. In 
the following, we focus on the properties of these four sources 
and primarily refer to them by their IGR names (even though 
Chandra names are also given to the sources). 

3. CHANDRA POSITIONS AND ASSOCIATIONS WITH KNOWN 
SOURCES 

Table|2]provides the Chandra positions for the four sources. 
The positions come from the wavdetect analysis described 
above, and the position uncertainties are 0".6. This value is 
the nominal Chandra pointing accuracy, and the statistical er- 
rors calculated by wavdetect are negligible in comparison. 
Here, we use the SIMBAD database to determine if the Chan- 
dra sources are associated with previously known sources. 

For IGR J16195^945, the Chandra source CXOU 
J161932.2-494430 lies I'.l from the best ASCA position 
for AX 1161929^945, which is just outside the 1' ASCA 
error circle (see Figure 0. While .Sugizaki et al.. (200 11) 
quote ASCA position uncertainties of 1', a careful analysis by 
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lUedaetd] shows that the 90% confidence error radii 

foTASCA sources are in the range 0'.6-0'.8, depending on the 
detection significance. Thus, the 1' error circle can be con- 
sidered to be sUghtly larger than a 90% confidence error cir- 
cle. Still, CXOU J161932.2-494430 and AX J161929-4945 
are close enough to consider the association likely. An asso- 
ciation between IGR J16195-4945, AX J161929-4945, and 
the B1/B2 la supergiant HD 146628 has also been suggested 
(jSidoli et al. 2005; Tomsick et al. 2004; Sugizaki et al. 200 1|. 
However, the Chandra position is 1'.2 away from HD 146628, 
clearly ruling out this association. There are no sources in the 
SIMBAD database within 1' of the Chandra position. 

There is no evidence that IGR J 16207-5 129 was previ- 
ously detected in the X-ray band. While the source is bright 
enough to be seen in the ASCA Galactic Plane Scans, it lies 
just outside of the region that was observed at a Galactic 
latitude of = -1.05°. A possible association between IGR 
JI6207-5129a nd the Al IVe star HD 146803 was suggested 
JTomsick et al.i r2004: Masetti et al. 2006b); however, this is 
ruled out by the fact that the positions for HD 146803 and 
CXOU J162046.2-513006 are not compatible. There are no 
sources in the SIMBAD database within 1' of the Chandra 
position. 

For IGR J16167^957 and IGR J 17 195-4100, the Chan- 
dra positions confirm associations with IRXS J161637.2- 
495847 and IRXS J171935.6-410054, respectively (see Fig- 
ure A SIMBAD search at the Chandra position of CXOU 
J161637.7-495844 reveals only IGR J16167-4957 and the 
ROSAT source within 1'. For IGR J 17 195-4100, the only 
other SIMBAD source within 1 ' of CXOU J 1 7 1 935 .8^ 1 0053 
is USNO-Bl.O 0489-0511283. This USNO source is 0".70 
from the Chandra position and must be considered as a possi- 
ble optical counterpart. We investigate this possibility further 
in §6. 

4. CHANDRA ENERGY SPECTRA AND LIGHT CURVES 

We produced 0.3-10 keV energy spectra for each of the four 
bright Chandra sources using the CIAO tool dmextract to 
produce the spectrum and the tools mkacisrmf and mkarf 
to produce the files needed to characterizing the instrument 
response. We used a circular source extraction region with a 
radius of 2". 5 and an annular background region centered on 
the source with inner and outer radii of 10" and 60", respec- 
tively. Given that two of the sources are close to the gaps be- 
tween the ACIS chips, we produced exposure maps to deter- 
mine the effective exposure time for each source. The sources 
IGR J 16207-5 129 and IGR J 17 195-4100 (ObsIDs 5472 and 
5474) received the full exposure time, but IGR J 16 195-4945 
(ObsID 5471) and IGR J 16 167-4957 (ObsID 5473) received 
70% and 67% of the full exposure time, respectively. The 
CIAO tool mkarf automatically accounts for the reduced in- 
strument response near chip boundaries. 

We fitted the spectra using the software package XSPEC 
11.3.2o. Originally, we rebinned the spectra and fitted them 
using x^-minimization. After rebinning, the spectrum for the 
lowest count rate source (IGR J16195-4945) had 9 bins with 
an average of 20 counts per bin, while we used 12-bin spec- 
tra for the other 3 sources. We fitted the spectra using a model 
consisting of an absorbed power-law, and, based on the results 
for other IGR sources, it would not be surprising if the col- 
umn density has intrinsic and interstellar contributions (e.g., 
[Walter et al. 2003; Revnivtsev ^^03 ). The absorbed power- 
law fit is acceptable (x^/v = 4.0/6) for IGR J 16 195^945, but 
we obtained very poor fits for the other three sources with val- 
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Fig. 2. — Chandra energy spectra and residuals for the four IGR sources fit- 
ted with an absorbed power-law model. For IGR J16207-5129, IGR J16167- 
4957, and IGR 117195^100, a pile-up correction is included in the model. 



ues of between 46 and 83 for 9 degrees of freedom (dof). 
Large positive residuals above ^^6 keV are the main reason 
for th e poor fits. This suggests that photon pile-up dDavisI 
1200 Ih is affecting these three spectra, and, indeed, one expects 
pile-up to significantly impact the spectral shape for sources 
with count rates as high as we observe for the three brighter 
sources (0.13-0.18 s~'). We refitted the three affected spectra 
after including the XSPEC pileup model (D aviSt2001 ). and 
the quality of the fits improves enormously with values of 
3.7, 8.7, and 9.2 (8 dof) for IGR J16207-5129, IGR J16167- 
4957, and IGR 117195^100, respectively. Comptonization 
models (e.g., comptt in XSPEC) approximate a power-law 
shape in the 0.3-10 keV band and provide fits of similar qual- 
ity to the power-law fits. While a thermal blackbody model 
can adequately describe the spectrum of the fainter source, 
such a model provides very poor fits to the spectra of the three 
brighter sources. Even if the pile-up correction is included, an 
absorbed blackbody model gives x^ values between 32 and 
98 for 8 dof. We conclude that these three spectra are non- 
thermal, and we focus on the power-law model below. 

Another signature of photon pile-up is the flattening of the 
point spread function (PSF), and this causes the three brighter 
sources to appear extended when their radial profiles (i.e., sur- 
face brightness vs. angle from the source position) are com- 
pared to the profiles of point sources that are unaffected by 
pile-up. We studied the radial profiles of the 4 IGR sources 
as a check of the pile-up correction, and we conclude that the 
PSF distortions are well-explained by pile-up only and that 
the pile-up corrections we apply in performing the spectral 
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Fig. 3. — Error contours for the parameters Nh and F from the fits to the 
Chandra spectra. The inner contour is 68% confidence (AC = 2.3) and the 
outer contour is 90% confidence (AC = 4.6). 



fits are correct. 

Figure 12 shows the 0.3-10 keV spectra and residuals for 
the absorbed power-law fits (with the pile-up correction), and 
Table provides the values of the spectral parameters. To 
obtain these values, we refitted the spectra by minimizing the 
Cash "C"-statistic (Cash 1979), which is an improvement over 
X^-minimization for spectra with small numbers of counts as 
the Cash analysis does not require discarding information by 
rebinning the spectra. While none of the sources show the ex- 
tremely high column densities that have been seen for some of 
the IGR sources, the measured values of A^h for IGR J 16 195- 
4945 and IGR J16207-5129 are 7!^ x lO^^ ^nd 3.7!} ^ x lO^^ 
cm"^ (90% confidence errors), placing them somewhat above 
their respective Galactic column densities o f 2.2 x 10^^ and 
1.7 X 10^2 cm-2 dPickev & Lo ckman"199ff). IGR J16167- 
4957 and IGR J17195-4100 have values of A^h that ai-e sig- 
nificantly lower than their Galactic column densities. 

All the sources are intrinsically hard in the 0.3-10 keV band 
based on the measured values of the power-law photon index. 
As shown in Table |3 the best fit values of F range from 0.5 
to 1.1, and, in all cases, the spectra are constrained so that 
F < 1.6 based on the 90% confidence (AC = 2.7) error bars. 
While the errors on the individual parameters indicate very 
hard spectra, we also calculated error contours to account for 
possible correlations between A^h and F. The contour plots 
in Figure 13 show the 68% (AC = 2.3) and 90% (AC = 4.6) 
confidence error contours, and it is clear that these param- 
eters are significantly correlated for all 4 spectra. While this 
analysis shows that the error regions are somewhat larger than 




5000 



Fig. 4. — 0.3-10 keV Chandra light curves with 64 s time bins for two of 
the IGR sources. IGR J16207-5129 is clearly variable, and there is evidence 
that IGR 117195^100 is also variable (see text). The light curves for the 
other two IGR sources are significantly affected by dithering over the gap 
between ACIS chips and are not shown. 



indicated by the individual parameter values given in Table|3l 
the spectra are still constrained to be quite hard, with F < 1 .8 
(90% confidence), and it is still the case that IGR J 16 195- 
4945 and IGR J 16207-5 129 have column densities that ai-e 
higher than the Galactic values while A^h is significantly less 
than the Galactic values for the other two sources. 

We produced 0.3-10 keV light curves with 64 s time bins 
for all 4 sources. The Hght curves for IGR J 16 195-4945 and 
IGR J 16 167^957 are strongly affected by the dithering over 
gaps between the ACIS chips. In both cases, the count rates 
drop to zero several times during the observation. In all, about 
30% of the time bins are affected, and we do not study the 
light curves for these 2 sources further. On the other hand, 
there are no problems with the light curves for IGR J 16207- 
5129 and IGR J 17 195-4100, and these ai-e shown in FigureH 
The clearest variability is seen for IGR J16207-5129. In the 
middle of the observation, the count rate drops to about a 
third of its average level, remains there for ^^500 s, and then 
recovers. We formally tested for variability by fitting each 
light curve with a constant and obtained x' /v = 163/78 and 
xVi' = 100/72 for IGR J 16207-5 129 and IGR J 17 195-4100, 
respectively. Thus, it is likely that both sources exhibit some 
variability. 

5. INTEGRAL ENERGY SPECTRA 

We also produced >22 keV spectra from the IBIS/ISGRI 
instrument on INTEGRAL. The observations used to produce 
the spectra include all public data available from INTEGRAL 
revolutions 46 to 235 (2003 February 27 to 2004 September 
15). We made sky images and mosaics with the Off-line- 
Analysis Software version 5.1 in 10 energy bands between 
22 keV and 250 keV, and we derived average source spec- 
tra from the mosaic images. After making vignetting, off- 
axis, and dead-time corrections, the effective exposure time 
is 495 ks for IGR J 17 195-4100 and -165 ks per source for 
the other 3 sources (see Table|3i. Despite the relatively long 
exposure times, these sources are rather faint for INTEGRAL 
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with 20-50 keV ISGRI count rates between 0.18 ±0.03 and 
0.35 ±0.03 s"' (compared to, e.g., a 20-60 keV count rate 
of 160 s~' for the Crab Pulsar and Nebula). Thus, the spec- 
tra of the IGR sources consist of a relatively small number 
of energy bins. For IGR J16195-4945, IGR J16207-5129, 
and IGR 117195^100, we are able to obtain 4 energy bins in 
which the source is detected at a significance of at least 4-ct, 
and, for IGR J 16 167^957, the source is only detected in 2 
energy bins. 

We fitted the ISGRl spectra with a power-law model, and 
the results of these fits are given in Table|3 For IGR J16195- 
4945, IGR J16207-5129, and IGR J17195^100, acceptable 
fits are obtained with the power-law model with values of 
r equal to 1.7 ±0.8, 1.9 ±0.5, and 2.8 ±0.8, respectively. 
While, for IGR 116195^945, the 90% confidence error re- 
gion for r overlaps with the error region determined using the 
Chandra spectrum, the spectra for IGR J 1 6207-5 1 29 and IGR 
117195^100 are significantly steeper in the 20-50 keV band 
than in the 0.3-10 keV Chandra band. Although the spectrum 
for IGR J 16 167^947 only has 2 energy bins, it is still possi- 
ble to calculate the value for F, and we obtain 4.3;!;] 4, which 
is significantly steeper than the Chandra measurement. 

The fact that the spectra are steeper in the ISGRI 20-50 
keV band when compared to the Chandra 0.3-10 keV band 
suggests that there is a high energy cut off. To investigate 
the cut off and to compare the flux levels measured by the 
two instruments, we performed combined fits for the Chandra 
plus ISGRI spectra. As there is evidence that these sources 
are variable, and the Chandra and ISGRI spectra were not 
taken simultaneously, it is essential to leave the overall nor- 
malization between the two instruments as a free parameter 
Leaving all parameters free, we fitted the spectrum for each 
source with and without an exponential cut off (highecut 
in XSPEC). The cut off is detected at the highest level of sig- 
nificance for IGR J 1 6207-5 1 29 with an improvement in the fit 
from x^/^= 19.6/1 1 to 5.3/9, corresponding to an F-test sig- 
nificance of 99.7%. The cut off is also marginally statistically 
significant for IGR J16167-4957, with a 98.4% chance that 
the cut off is required. The spectra for IGR J 16 195-4945 and 
IGR J 17 195^100 aflow for the possibihty that there is a cut 
off in the ISGRI energy range, but the cut off is not required. 

While we allowed the A^h, the power-law parameters, a 
in the pile-up model (for the three brighter sources), the cut 
off parameters (E^ut and iifoid), and the overall normaliza- 
tion between the two instruments to be free for the fits de- 
scribed in the previous paragraph, with so many free parame- 
ters, the parameter values are very poorly constrained. Thus, 
we fixed the A^h, power-law, and pile-up parameters to the 
values measured by Chandra, and left just the normaliza- 
tion between instruments and the exponential folding energy 
(£^foid) as free parameters in the fit. We fixed the energy at 
which the exponential cut off begins (ficut) to 10 keV, just 
above the Chandra bandpass. However, it should be noted 
that, while the Chandra spectra do not show any evidence for 
a spectral cut off, the quality of the spectra and the distor- 
tion of the spectra due to pile-up do not allow us to rule out 
the possibility that the cut off starts at a lower energy. For 
IGR J16195-4945, IGR J16207-5129, IGR J16167-4957, 
and IGR J 17 195-4100, respectively, we obtained ISGRI rel- 
ative normalizations of l.l![;;^, 1.4+g|, and 0.9!j]:^ 
as well as the following values for Em^: 32'!i\l^, 23;!;^-', 9^1, 
and 20;!;^^ keV (90% confidence errors). These values indi- 
cate that IGR J 16207-5 129 was relatively bright during the 



Chandra observation compared to the brightness detected by 
ISGRI in 2003-2004, while the other sources were consistent 
with a relative normalization of 1.0. Also, it is notable that 
the range of values we measure for Efou, 9-32 keV, are simi- 
lar to the values that have been measured for HMXB pulsars 
(iCoburn et al..2002i) . 

6. IR AND OPTICAL IDENTIFICATIONS 

With the 0".6 Chandra positions for the 4 IGR sources, 
we are able to search for counterparts to these sources at 
other wavelengths. We searched the following catalogs: 
United States Naval Observatory (USNO-Bl.O and USNO- 
A2.0); Deep Near Infrared Survey of the Southern Sky (DE- 
NIS); 2 Micron All-Sky Survey (2MASS); and Spitzers 
Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE). We also performed optical and IR photometry 
at ESO's New Technology Telescope, and the NTT observa- 
tions and results are described below. 

The USNO catalogs are from optical surveys that cover 
most of the sky. For the USNO-Bl.O, the position accura- 
cies are 0".2 and the photometry is good to ±0.3 magnitudes 
in B, R, and /. USNO-A2.0 provides 0".25 positions and B 
and R measurements that are generally accurate to ±0.5 mag- 
nitudes but may be somewhat worse that this in the southern 
hemisphere. DENIS^ is an optical and IR survey, concentrat- 
ing on the southern hemisphere. The position accuracies are 
1", and accurate photometry is available in the /, /, and K, 
bands. 2MASS** is an all-sky IR survey. The catalog provides 
very good positions (0".2 accuracy) and accurate photometry 
in the J, H, and bands. GLIMPSE uses the Infrared Ar- 
ray Camera (IRAC) on Spitzer to obtain short exposures of 
the Galactic plane, covering the Galactic longitude range be- 
tween 10° and 65° on both sides of the Galactic center and 
Galactic latitudes within 1 ° of the plane. Images are obtained 
in four bands centered at wavelengths of 3.6, 4.5, 5.8, and 8.0 
fim. The GLIMPSE team has compiled an on-line catalog' 
with lists of detected sources, including fluxes and positions 
accurate to 0".3. 

On 2004 July 10 between UT and 1 hr, we obtained op- 
tical photometry in B, V, R, and / bands of the fields of IGR 
J 16 195-4945 and IGR J 16207-5 129 with the spectro-imager 
EMMI installed on the NTT. We used the large field imaging 
of EMMI's detector, giving, after 2x2 rebinning, a pixel size 
of 0".333 and a field of view of 9'.9 x 9'.0. We used an inte- 
gration time of 1 s per exposure. We observed 5 pho tometric 
standard stars from the optical standard star catalog ol'LandolJ 
( 1992): PG 1633H-099, PG 1633H-099A, PG 1633H-099B, PG 
1633H-099C, and PG 1633H-099D. We also obtained IR pho- 
tometry in J, H, and K, bands of IGR J 16207-5 129 (on 2004 
July 8 at UT 4.67 hr), of IGR J 16 167-4957 (on July 1 1 at UT 
5.67 hr), and of IGR J 17 195-4100 (on July 11 at UT6.83hr) 
with the spectro-imager SofI installed on the NTT. We used 
the large field imaging of Sofl's detector, giving a pixel size 
of 0".288 and a field of view of 4'.94 x 4'.94. These photo- 
metric observations were obtained by repeating a set of im- 
ages for each filter with 9 different 30" offset positions, in- 
cluding the targets, with an integration time of 60 s for each 
exposure, following the standard jitter procedure allowing for 
clean subtraction of the blank sky emission in IR. On several 
occasions, we observed 3 photometric standard stars from the 

^ See http://cdsweb.u-strasbg.fr/denis.html 

See http://www.ipac. caltech.edu/2 mass/reIeases/second/doc/| 
' See http://www.astro.wisc.edu/sirt^] 
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Fig. 5.— NTT images for(a) IGR 116195^945 (/-band), (b) IGR J16207- 
5129 (/f,-band), (c) IGR J16167^957 (K,-hand), and (d) IGR 117195^100 
(^Ts-band). The lengths of the N and E an'ows are 2". From smallest to 
largest, the error circles are: 2MASS (thin white, 0".2 radius error circle), 
USNO-Bl.O (thick grey, 0".2), USNO-A2.0 (thick white, 0".25), GLIMPSE 
(thin grey, 0".3), Chandra (thin white, 0".6), DENIS (thin grey, l"). 



faint IR standard star catalog of lPersson et alJJl998h : sj9157, 
sj9172, and sj9181. We used the Image Reduction and Anal- 
ysis Facility (IRAF) suite to perform data reduction, carrying 
out standard procedures of optical and NIR image reduction, 
including flat-fielding and NIR sky subtraction. We carried 
out aperture photometry, and we then transformed instrumen- 
tal magnitudes into apparent magnitudes. The targets were 
close to the zenith during the observations. 

Table |5] provides the results of our search of the catalogs, 
and also the results of our NTT observations. For each of 
the 5 catalogs, we list the closest source to the Chandra posi- 
tion and the measured magnitudes. In Figure |5] we show an 
/-band NTT image for IGR 116195^945, and K.-hand im- 
ages for the other three sources. The Chandra positions and 
the positions obtained from the catalogs are plotted on these 
images. 

For IGR J16195-4945, the Chandra position lies on an 
apparently stellar IR source (based on an inspection of the 
2MASS images) with a magnitude of 7 = 13.57 ± 0.03. The 
2MASS and GLIMPSE positions lie inside the Chandra error 
circle, indicating a good association between the X-ray and IR 
sources (see Figure (5^). However, the USNO-A2.0 position 
lies on the edge of the Chandra error circle, and the USNO- 
Bl.O position lies slightly outside. This may indicate the pres- 
ence of a second optical source blended with the counterpart. 
Due to the faintness of the source (or sources) it is not clear 
from the NTT /-band image whether there are two blended 
sources or not. However, we present analysis of the spectral 
energy distributions (SEDs) below, and the IGR J 16 195-4945 
SED (see Figure|6t) supports the possibility that there are two 
blended sources. 

For IGR J 16207-5 129, the Chandra position lies on a 
bright IR source with J = 10.44 ± 0.02. This source is present 
in all 5 catalogs, and all the positions are consistent with the 
Chandra source except for the USNO-Bl.O position, which 
lies just outside the Chandra error circle (see Figure [Sj)). 



While this could indicate the presence of an interloper, the 
SED does not suggest any contamination by other sources (see 
Figure I^J)). 

For IGR J16167^957, the Chandra position lies on a 
source with J = 14.86 ± 0.06. The source appears in the 
GLIMPSE, 2MASS, and USNO catalogs, and all positions 
lie within the Chandra error circle (see Figure|5j;). Although 
the 2MASS images suggest that the source is extended to 
the South-East, the NTT /Tj-band image shows that there is 
actually a blend of at least 3 sources. The IGR J 16 167- 
4957 Chandra position is consistent with the brightest of the 
blended IR sources, but there is some suggestion in the NTT 
image that the brightest "source" may also be a blend. 

For IGR J 17 195^100, the Chandra position is near two 
sources that are clearly blended in the 2MASS images. The 
position of the 2MASS source is 2". 6 away from the Chan- 
dra position, indicating that the bright IR source is not IGR 
117195^100. With /-band and IR coverage, both sources ap- 
pear in the DENIS survey, and the optically-brighter, North- 
West source is the one that is consistent with the Chandra 
position. Our results confirm the association between IGR 
J 1 7 1 95-4 1 00 and USNO-B 1 . 0489-05 11283. 

7. SPECTRAL ENERGY DISTRIBUTIONS 

Figure|6lshows the spectral energy distributions (SEDs) for 
the four sources. These include the optical and IR measure- 
ments given in Table |5] the spectrum in the Chandra band- 
pass, and INTEGRAL measurements. It should be noted that 
the optical and IR fluxes are not dereddened, and absorption 
is not removed from the Chandra spectra. The SED measure- 
ments are spread over a period of years, so discontinuities in 
the SED may be due to source variability. For the IR and 
optical parts of the SED, a check on variability is provided 
by the fact that we have multiple measurements for many 
of the photometric bands. A comparison shows that the IR 
fluxes are quite stable for IGR 116195^945, IGR J16207- 
5129, and IGR J16167-4957. For IGR J 16 195-4945, the 1-cr 
error regions overlap for the J- and /STj-band measurements 
from DENIS and 2MASS. For IGR J16207-5129, the DE- 
NIS, 2MASS, and NTT J- and K.-hund measurements only 
range from 10.38-10.54 and 9.13-9.18,respectively. For IGR 
116167^957, the 2MASS and NTT measurements are con- 
sistent at the 1-0- level for J, H, and K,. For IGR J 16207- 
5129, the optical flux also appears to be rather stable, but for 
IGR J 16 195^945 and IGR J 16 167-4957, the agreement be- 
tween the various optical measurements is not as good (see 
Figure|6j. For IGR J17195-4100, the agreement in the opti- 
cal bands is relatively good, but this source has the sparsest 
optical/IR measurements. 

For IGR 116195^945, we simply show in Figure [6^ the 
data and the best fit model for Chandra; however, for the other 
three sources, the spectrum is significantly distorted by pile- 
up, and we show the range of flux measurements based on 
the error ranges for the model parameters. For IGR J16207- 
5129, the flux range shown reflects the 90% confidence error 
range on L, which is 0.0-1.1. For IGR J 16 167-4957 and IGR 
J17195-4100, the T ranges are 0.7-1.6 and 0.8-1.4, respec- 
tively. The INTEGRAL 20-50 keV measurements and 50-250 
ke V upper limits are shown. As described in §5, the fits shown 
include a high energy cut off and a free normalization between 
Chandra and ISGRI. IGR J 16207-5 129 is the only source for 
which the normalization for ISGRI relative to Chandra is re- 
quired to be less than 1 .0, and this is apparent in Figure|6l 

To characterize the shape of the optical and IR contin- 
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Fig. 6. — Spectral energy distributions (SEDs) for tlie four IGR sources. 
The IR and optical fluxes come from the measurements given in Tabled and 
this part of the spectrum has been fitted with a blackbody model modified by 
the effects of extinction. The dotted line indicates the best fit (see text). Chan- 
dra and INTEGRAL measurements are also shown. For IGR J16207-5129, 
IGR J16167^957, and IGR J17195^100, our best estimate of the specQ^m 
in the Chandra energy range is the model (solid line). For the optical/IR data, 
the different symbols are: GLIMPSE (filled squares), 2MASS (open squares), 
DENIS (filled triangles), USNO-Bl.O (filled diamonds). USNO-A2.0 (open 
diamonds), NTT (filled circles). 



uum, we fitted each of the SEDs with four models: a sin- 
gle blackbody, a single power-law, the sum of 2 blackbod- 
ies, and a blackbody plus a power-law. Presumably, black- 
body components could be indicative of thermal components 
from a star, but we consider the power-law to be a phe- 
nomenological component. The models also account for 
interstellar extinction using the analytical approximation of 
ICardelh. Clavton & MathisI (11989*). This adds a single free 
parameter. Ay, to the single component models, while, for 
the two-component models, we allow the components to have 
different extinctions. We used x^-minimization when fitting 
the SEDs, and, for IGR J 16 195-4945 and IGR J 16207-5 129, 
all of the models resulted in formally very poor fits with the 
best models leading to x^/j^ =128/14 and x^/^ = 684/16 for 
the two sources, respectively. Such poor fits may indicate in- 
adequacies of our models, but they are also at least in part 
due to systematic errors, which may be related to the vari- 
ability discussed above or possibly due to, e.g., uncertainties 
in magnitude/flux conversion. While it is clear that system- 
atic errors must be included for these 2 sources, it is not clear 
which data points are the most suspect. Thus, we included 
systematic errors by simply multiplying the existing errors on 
all of the optical and IR data points by a constant factor. We 
adjusted the constant factor until we achieved a reduced-x^ 



of 1.0 for at least one of the models. For IGR 116195^945 
and IGR J16207-5129, the constant factors are 2.2 and 6.5, 
respectively. For the other 2 sources, IGR J 16 167^957 and 
IGR J 17 195^100, we did not need to include any systematic 
errors. The SED fit results are presented in Table|6l and here, 
we describe the results for each source in turn. 

For IGR J 16 195-4945, neither of the single component 
models provide an acceptable fit to the optical/IR SED, and 
a two-component model is required at high significance. The 
two blackbody and blackbody plus power-law models provide 
acceptable fits, and both imply that the IR flux is dominated 
by a high temperature star (>9400 K) with high extinction 
(Ay ^ 18). Also, both models imply that the component that 
dominates the optical flux has a significantly lower extinction 
(Ay < 7.1). This SED and the shift between the optical and 
IR source positions reported in §6 provide evidence that the 
optical component comes from an interloping star, while the 
IR component likely comes from IGR J 16 195-4945. 

Our highest quality optical/IR SED is the one for IGR 
J16207-5129. For this SED, two-component models do not 
provide any improvement over the single component mod- 
els, and the blackbody and power-law models provide fits 
of roughly the same quality. The power-law index is a = 
-2.29;!;q 2g, which is nearly consistent with the value of a = -2 
expected for a blackbody. The fits strongly suggest that the 
spectrum is thermal as would be expected for stellar emission. 
From the blackbody fit, the value obtained for the extinction. 
Ay = 10.8;!;Q g is very close to the value of Ay = 9.5 obtained 
by converting the Galactic column density of A^h = 1 .7 x 10^^ 
cm~^ to optical extinction using the lPredehl & Schmittl(IT995h 
relationship. 

A blackbody model does not provide a good fit to the IGR 
J 16 167-4957 and IGR J 17 195^100 SEDs with values of 
X^/i^ of 54/10 and 26/6, respectively. However, both SEDs 
are adequately described by a power-law model with indices 
of a = -0.36!l;|j [g and a = 0.17;!;q 47, respectively. The power- 
law fits for both sources imply low extinction, which is con- 
sistent with the values of A^h measured via fits to the X-ray 
spectrum. Although we fitted the SEDs for these two sources 
with the two component models, the improvements in x' do 
not justify the addition of another component, and the param- 
eters for the two-component fits are very poorly constrained. 

8. DISCUSSION 

In the following, we discuss the detailed constraints on 
the nature of each source, considering the information we 
have obtained from Chandra, INTEGRAL, and the optical/IR 
SEDs. 

8.1. IGR 116195^945 

Of our 4 sources, IGR J16195-4945 is the most similar to 
the class of obscured IGR sources. At A^h = (7^3) x 10^^ cm"^, 
the X-ray measured column density is higher than the Galac- 
tic value of 2.2 X 10^^ cm"^, suggesting the possibility that the 
X-ray source is intrinsically absorbed. Due to the possibility 
of intrinsic absorption, we cannot use the X-ray-measured A^h 
to constrain the distance to the source; however, because the 
optical and IR emission from the system is likely to predomi- 
nantly come from regions that are not as affected by intrinsic 
absorption, such as the companion star, we expect that the op- 
tical extinction is dominated by interstellar material. Thus, 
the fact that Ay = 17.5;!;2 3 ^'^^ '■^^ extinction on the blackbody 
component is as large or larger than the value obtained by 
converting the Galactic A^h to optical extinction (Ay = 12.3), 
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Fig. 7.— Error contours for IGR 116195^945 (a) and IGR J16207-5129 
(b) for the blackbody parameters in the opticalAR SEDs. Contours are shown 
for the stellar radius (R) and temperature (T), assuming various distances. 
The solid contours are 68% confidence (Ax^ = 2.3), and the dashed contours 
are 90% con fidence (Ay^ = 4.6). Radii and temperatures for various spectral 
types from lde Jager & Nieuwenhuiizea Iil987) are overplotted, allowing us 
to constrain the likely spectral types of the companions in these systems. 



is likely to indicate that we are looking at the source through 
much of the Galaxy. The distance to the Norma-Cygnus arm 
along this line of sight is ^5 kpc ( Russeil 2003 ), and, although 
we cannot formally constrain the distance to the source, we 
assume a fiducial value of 5 kpc. From the unabsorbed fluxes 
reported in Table |3] and 0] the implied X-ray luminosities are 
1.4 X lO^'* (d/5 kpc)2 (0.3-10 keV) ergs s"' and 5.8 x 10^"* 
(d/5 kpc)2 ergs s"' (20-50 keV). 

The IR portion of the SED is consistent with a single- 
temperature blackbody, suggesting that it is more likely that 
the IR emission comes from a stellar companion than from, 
e.g., an accretion disk. The lack of variability in the IR emis- 
sion discussed in §7 also provides evidence that the emis- 
sion comes from a star. The fit to the SED gives a lower 
limit on the blackbody temperature of T > 9400 K, and 
this indicates a spectral type earlier than A1-A2 regardless 
of whether the star is a main sequence star or a supergiant 
(|de Jaaer & Nieuwenhuiizen 1987; Cox 2000). The question 
of whether the companion is a main sequence star or a super- 
giant depends sensitively on the distance to the source. In the 
fit to the SED, the parameters T and R/d have a significant 
level of degeneracy, and, in FigureQi, we show 68% and 90% 
confidence error contours for T and R for assumed distances 
of 5 kpc, 10 kpc, and 15 kpc along with temperatures and radii 
for stars with various spectral types. If the distance is ^3-9 
kpc, the companion is a main sequence B or O-type star If 
flie distance is ^^9-15 kpc, a large range of supergiant A, B, 
and O-type spectral types are possible. 

8.2. IGR J16207-5129 

Like IGR 116195^945, the IGR J 16207-5 129 A^h is some- 
what (about a factor of 2) higher than the Galactic value, al- 
lowing for, but not requiring, some intrinsic absorption. From 



the optical/IR SED, Ay = 10.8;!;o g, ^^'^ ^^^^ value suggests a 
relatively large distance for this source. Thus, it is reasonable 
to take 5 kpc, the distance to the Norma-Cygnus arm along the 
line of sight, as a fiducial distance. The implied 0.3-10 keV 
luminosity is 1.3 x lO-'^ {d/5 kpc)^ ergs s"', about an order of 
magnitude higher than IGR 116195^945, and the 20-50 keV 
luminosity is 8.3 x lO-''* {d/5 kpc)^ ergs s"^ 

In this case, the optical/IR SED consists of a single- 
temperature blackbody, and it is very likely that this emission 
comes from the stellar companion. The lower limit on the 
stellar temperature is > 18,000 K, indicating that this system 
harbors a very hot, and likely massive, star Figure Q shows 
the error contours for the parameters T and R assuming dis- 
tances of 5 and 10 kpc. At ^^3-5 kpc, the companion is a main 
sequence B or O-type star At ^5-9 kpc, various supergiant 
types are possible. At 10 kpc, the range of possible temper- 
atures and radii are not consistent with any known spectral 
type, indicating that this system is closer than 10 kpc. 

8.3. IGR J16167^957 

Although this source is close to the Galactic plane at = 
+0.50° and is within 0.5° of IGR J16195-4945 and 1.7° of 
IGR J 16207-5 129, its X-ray spectrum indicates a column 
density of A^H = (5!|)x 10^' cm , which is an order of mag- 
nitude lower than the value for IGR 116195^945 and is also 
lower than the Galactic column density of 2.2 x 10^^ cm"^ 
along this line of sight. The measured A^h corresponds to an 
extinction of Ay = 2.8^} j, and, although the optical/IR SED 
is fit with a phenomenological power-law model, the implied 
extinction is even lower. Ay = 13^0 5. While the inhomoge- 
neous distribution of gas and dust in the Galaxy make it diffi- 
cult to use extinction measurements to determine source dis- 
tances in any rigor ous way , here, we make an estimate using 
the values from the lDiplas & Savage ( 1994) study, which in- 
cludes UV extinction measurements of nearly 400 early-type 
stars through various lines of sight through the Galaxy and 
with distances as large as 10 kpc. In that study, the average 
extinction per kpc is Ay /d = 0.45 magnitudes/kpc. For IGR 
116167^957, the extinction inferred from the X-ray spec- 
trum allows the distance to be anywhere from 2 to 10 kpc, 
but the optical/IR extinction implies a distance of 2.9;!;['[ kpc. 
This estimate is marginally consistent with the source being 
in the Norma-Cygnus arm at ~5 kpc, but the line of sight also 
crosses the Scutum-Crux arm at ~3 kpc. Given the lower ex- 
tinction for this source compared to IGR 116195^945 and 
IGR J 16207-5 129, which ai-e likely in the Norma-Cygnus 
arm, we take 3 kpc as a fiducial distance for IGR J 16 167- 
4957, and the X-ray luminosities are 4.7 x lO"*"* {d/3 kpc)^ 
ergs S-' (0.3-10 keV) and 1.7 x 10^"* {d/3 kpc)^ ergs s"' (20- 
50 keV). 

The interpretation for the shape of the optical/lR SED is not 
immediately clear, but the fact that it is not consistent with 
a single temperature blackbody suggests that the emission is 
not simply from a companion star. The power-law index of 
-0.36;!;Q [g is considerably flatter than a blackbody, and this 
could indicate a multi-temperature blackbody as might be ex- 
pected from an accretion disk. Regardless of the dominant 
contributor to the optical/IR emission, we can take the mea- 
sured magnitudes as an upper limit on the contribution from 
a putative optical companion. If we assume a distance and 
an extinction at the upper ends of the derived ranges, 4.2 kpc 
and Av = 2.8, respectively, and interpolate between the mea- 
sured B and R magnitudes to infer a V-band magnitude of 
16.5, the limit on the absolute magnitude is My > 0.6. Unless 
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the distance is very much larger than 4.2 kpc, this rules out 
a supergiant companion, and the companion's spectral type is 
later than AOV. Thus, IGR J 16 167-4957 is not an HMXB. 

8.4. IGRJ17195^100 

At I = 346.98° and b = -2.14°, this source is -14° from 
the other 3 sources, and is slightly off the Galactic plane. The 
measured upper limit on the column density of A^h < 2 x 10^' 
cm~^ is considerably lower than the Galactic value of 7.7 x 
10^' cm~^, indicating that the distance to the source may be 
relatively low. The A^h from the X-ray spectrum implies Ay < 
1.18, and the optical/IR SED gives an up per limit of Ay < 
0.96. Using the average value ofAy/d from lDiplas & Savag d 
11994) indicates a distance upper limit in the range d <2.1- 
2.6 kpc. The upper limit on the distance of 2.6 kpc implies 
upper limits on the X-ray luminosities of <2.0 x lO-''* ergs s"' 
(0.3-10 keV) and <1.6 x lO^^gj-gs s"' (20-50 keV). 

Like IGR J16167^957, the physical interpretation for the 
power-law shape of the optical/IR SED is not clear. We can 
take the optical/IR magnitudes observed as upper limits on the 
contribution from a companion star. Interpolating between the 
B and /?-bands, we estimate that V = 14.8. This, along with 
the highest upper limits we measure for the extinction. Ay < 
1.18, and distance, d < 2.6 kpc, imply an upper limit on the 
absolute magnitude of My < 1.5. This rules out a supergiant 
companion or an OB main sequence star The earliest spectral 
type possible is '^A3V. 

9. SUMMARY OF THE NATURE OF THE 4 IGR SOURCES 

Table contains a summary of the results described in §8, 
including constraints on the source distance. X-ray luminosi- 
ties, and possible spectral types. Our basic conclusion is 
that IGR J 16 195-4945 and IGR J 16207-5 129 are very hkely 
HMXBs and are at distances consistent with the Norma- 
Cygnus arm (although the distances are not well enough con- 
strained to say definitively that they are in this spiral arm), 
while the other two IGR sources (IGR J 16 167-4957 and IGR 
117195^100) are nearer and are not HMXBs. Assuming dis- 
tances of 5 kpc for the 2 HMXBs, the 0.3-10 keV luminosities 
are lO^'*"^^ ergs s"', which is, not surprisingly, significantly 
lower than the brightest HMXB pulsars (e.g.. Vela X-1, Cen 
X-3), but the luminosities are not atypical when considering 
the full range of HMXB pulsar luminosities. 

Of the 4 IGR sources, only IGR J 16 195-4945 has a col- 
umn density that is consistent with the source having signifi- 
cant intrinsic absorption. While IGR J 1 6207-5 1 29 could con- 
ceivably also have some level of intrinsic absorption, com- 
paring its column density to the Galactic value suggests that 
the intrinsic A^h is not more than — 1 x 10^^ cm"^ while many 
IGR sources have A^h ^ 10^^"^'* cm~^. From the perspective 
of IGR sources and HMXB evolution jDean et a l.i .2005) . the 
very high column density IGR sources may be in an evolu- 
tionary stage during which the neutron star is spiraling to- 
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ward its companion, becoming embedded in the wind from 
the high-mass companion. In this picture, we may be see- 
ing IGR J 16207-5 129 early in this evolutionary phase so that 
the neutron star is still in a low density portion of the wind. 
Source-to-source differences in A^h could also be indicative of 
variations in the strength of the stellar wind over time. 

After submission of this work, two other reports provided 
further information about the HMXBs. Using INTEGRAL, 
Sgueraetal. (2006) found a 1.5 hour X-ray outburst from 
IGR J16195-4945, indicating that this source may belong to 
the class of supergiant fast X-ray transients. Also, for IGR 
J 16207-51 29, Masetti et al. (2006a) report on optical spectra 
of USNO-A2.0 0375-27093111 (i.e., the same star that we 
have determined to be the counterpart), and find that the spec- 
trum includes Ha in emission and is consistent with that of an 
HMXB. 

iMasetti et al. (2006a) also report on the optical spectra of 
the stars that we have (independently) determined to be the 
optical counterparts of IGR J 17 195-4100 and IGR J 16 167- 
4957. These stars show Hydrogen Balmer series emission 
lines as well as Hel and Hell emission lines. While this in- 
dicates that they are either Cataclysmic Variables (CVs) or 
LMXBs, the hard X-ray spectra would be very unusual for 
LMXBs (e.g., Muno et al. 2004). A CV interpretation is more 
likely, and iMasetti et alj(l2006al) also come to this conclusion. 
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TABLE 1 
Chandra OBSERVATIONS 



ObsID 


Target /" fo* Start Time 


Exposure Time (s) 


5471 


IGRJ16195^945 333.56 +0.34 2005 April 29, UT 17:25 


4,752 


5472 


IGR J 16207-5 129 332.46 -1.05 2005 June 25, UT 03:17 


5,109 


5473 


1GRJ16167^957 333.06 +0.50 2005 June 13, UT 19:14 


4,979 


5474 


IGR J17195^100 346.98 -2.14 2005 July 25, UT 12:57 


4,701 



"Galactic longitude in degrees. 
''Galactic latitude in degrees. 



TABLE 2 

Chandra COUNT RATES AND POSITIONS 



IGR Name 


CXOU Name 


ACIS Rate" 


R.A.* 


Decl.'' 




X-Ray Identification 


J16195^945 


J161932.2-494430 


0.039 


16''19°'32^20 


-49°44'30", 


.7 


AXJ161929^945(?) 


J 16207-5 129 


J162046.2-513006 


0.13 


16''20"' 46^26 


-51°30'06", 


.0 




J16167^957 


J161637.7^95844 


0.18 


16'M6'"37^74 


-49°58'44", 


.5 


IRXS J161637.2^95847 


J17195^100 


J171935. 8^10053 


0.15 




-41°00'53", 


.6 


IRXS J171935.6^10054 



"Count rate detected by the ACIS-I instrument in the 0.3-10 keV bandpass. 
''Right Ascension (equinox J2000). The radius of the eiTor circle is 0".6. 
'^Declination (equinox J2000). The radius of the error circle is 0".6. 



TABLE 3 
Chandra SPECTRAL RESULTS 



IGR Name 


(Xl022 cm-2) 


r 


Flux* 




Fit Statistic"' 


Galactic A'h'' 
(xlO^- cmr^) 


J16195^945 
J 16207-5 129 
J16167^957 
J17195^100 


n 

3.7!|;^ 

0.5ti 
0.08^;'i 


5+" '^ 

1 1+0.5 

1.1±0.3 


4 6+- ' 
^■"-0.8 

42!^ 

44+11 

25!| 


43+*"' 

"•"^-0.07 

64+<'"'' 
"■°^-o.io 


453 
667 
701 
743 


2.2 
1.7 
2.2 
0.77 



"Errors in this table are at the 90% confidence level (AC = 2.7). 
*Unabsorbed 0.3-10 keV flux in units of 10"'^ erg cm~ ^ s~* . 

^The grade migration parameter in the pile-up model IDavidlTOOll) . The probability that n events will be piled together but will still be retained after data 
filtering is o""' . 

''The Cash statistic for the best fit model. In e ach case, the spectra include 663 energy bins. 
''The column density through the Galaxy from lDickev & Lockmaniri990l) . 
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TABLE 4 
INTEGRAL Spectral Results 



ICR Name Exposure Time (ks) Count Rate" T Flux'' X'ii^ 

J16195^945 167 0.24 ±0.03 1.7 ±0.8 19 ±3 1.3/2 

J16207-5129 165 0.35 ±0.03 1.9±0.5 28±3 2.5/2 

116167^957 167 0.18 ±0.03 4.3!| J 16!^ 0/0 

117195^100 495 0.23±0.02 2.8±0.8 19±3 1.1/2 

"20-50 keV ISGRl count rate. 

*20-50 keV flux in units of 10"'^ erg cm-^ sT^ . 



TABLE 5 

Optical and IR magnitudes or fluxes 





IGR 116195^945 


IGRJ16207-5129 


IGR 116167^957 


IGR 117195^100 


USNO-Bl.O 


Name" 


0402-0529810 


0384-0560875 


0400-0527262 


0489-0511283 


Separation* 


1".03 


0".81 


0".53 


0".70 


B 


16.5 ±0.3 


19.7 ±0.3 


16.8±0.3 


15.1 ±0.3 


R 


15.3±0.3 


15.2±0.3 


15.8 ±0.3 


14.4 ±0.3 


I 


15.6±0.3 


13.0±0.3 


14.5 ±0.3 


14.3 ±0.3 


USNO-A2.0 


Name" 


0375-27014824 


0375-27093111 


0375-26829054 


0450-27095307 


Separation* 


0".72 


0".31 


0".45 


0".31 


B 


17.4 ±0.5 


18.9 ±0.5 


16.8 ±0.5 


15.2 ±0.5 


R 


16.2 ±0.5 


15.6±0.5 


16.2 ±0.5 


14.3 ±0.5 


DENIS 


Name" 


1161932.1^94430 


J162046.2-513006 




1171935.8^10053 


Separation* 


0".42 


0".29 




0".10 


/ 


15.38 ±0.05 


13.4±1.0 




14.81 ±0.04 


J 


13.55 ±0.08 


10.54 ±0.05 






Ks 


10.92 ±0.06 


9.17 ±0.06 






2MASS 


Name" 


116193220^944305 


J16204627-5130060 


J16163776-4958445 


J17193608-4100548 


Separation* 


0".14 


0".ll 


0".26 


2" .6 


J 


13.57 ±0.03 


10.44 ±0.02 


14.86 ±0.06 


C 


H 


11.96 ±0.03 


9.62 ±0.02 


14.28 ±0.09 


C 


K, 


11.00 ±0.02 


9.13 ±0.02 


13.76±0.10 


c 


GLIMPSE 


Name" 


G333.5571-i-00.3390 


G332.4590-01.0501 


G333.0560-h00.4962 




Separation* 


0".21 


0".07 


0".34 




3.6 /^m 


26.5 ±1.1 mJy 


100±3 mJy 


1.79±0.16mJy 




4.5 jua 


20.3 ±1.1 mJy 


61±3mJy 


1.40±0.16mJy 




5.8 /im 


14.3 ± 1.0 mJy 


46±2mJy 






8.0 //m 


8.5±0.5mJy 


28.4 ± 0.9 mJy 






New Technology Telescope 


B 


18.14 ±0.06 


19.8±0.1 






V 


17.22 ±0.05 


17.74 ±0.06 






R 


16.42 ±0.05 


15.38±0.03 






I 


15.54 ±0.03 


13.58 ±0.02 






J 




10.38 ±0.02 


15.0 ±0.1 


14.1 ±0.1 


H 




9.60 ±0.02 


14.4±0.1 


13.65 ±0.07 


K, 




9.18 ±0.04 


13.8±0.1 


13.2±0.1 



"Name of the nearest star in the catalog. 

*Separation between the Chandra position and the nearest star in the catalog. 
"^This 2MASS source is too far from the Chandra position to be associated. 
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TABLE 6 
IR/Optical SED Fits 





IGRJ16195^945 


IGRJ16207-5129 


IGRJ16167-4957 


IGR 117195^100 


Blackbody 


Av 


<0.16 


> 18000 
<5.53 


<1.16 


<2.84 


r(K) 

R/d (R0/kpc)" 


1944!3« 
10.3 ±0.6 


2790+35° 
1.84 ±0.21 




X /v 


269/16 


23/19 


54/10 


26/6 


Power-law 


Av 


7.1 ±0.6 
-0.73 ±0.13 
0.030 ±0.003 


1 1 6+"'* 

^•^'-0.28 

0.27 ±0.03 


1 3+0.6 

-0 36+0 " 

^•-'"-0.18 

0.00191 ±0.00015 


<0.96 
0.0039 ± 0.0006 




300/16 


19/19 


12/10 


9/6 


2 Blackbodies 












7l(K) 


>9400 








Ri/c;(R0/kpc) 


<4.36 










<4.89 








7-2 (K) 


>3800 








/?2/^/ (RQ/kpc) 


<1.13 










12/14 


23/17 


8/7 


7/3 


Blackbody+Power-law 




17.5!«? 








r(K) 


>9500 








Rjd (Ro/kpc) 


<4.25 








a° 


4.3«-8 
-1 7+i° 
OOlS+ooo^o 

12/14 


18/17 


8/7 


5/3 



"R/d if, the ratio of the radius of the spherical blackbody (in Rq) to the distance to the source (in kpc). 
''The power-law index, defined according to h], 3c . 

"^The power-law normalization, corresponding to the flux in Jy at a reference frequency of lO''' Hz. 



TABLE 7 
Summary of Results 



IGR Name 


d (kpc) rfflduciai (kpc) 


L (0.3-10 keV)" 


t (20-50 keV)* 


Spectral Type 


J16195^945 
J16207-5129 
J16167-4957 
J17195-4100 


~3-15 5 (Norma-Cygnus?) 
~3-10 5 (Norma-Cygnus?) 
1.8^.2 3 (Scutum-Crux?) 
<2.6 <2.6 


1.4 X 10^'' 
1.3 X 1035 
4.7 X 10^* 
<2.0 X lO^* 


5.8 X lO''' 
8.3 X 10'" 
1.7 X lO^" 
<1.6x 10'" 


OBV or supergiant 
OBV or supergiant 
later than AOV 
later than A3V 



"The X-ray luminosity in ergs s measured by Chandra at the fiducial distance. 
*The X-ray luminosity in ergs s"' measured by INTEGRAL at the fiducial distance. 



